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A new class of quantum channels called Grassmann channels is introduced and their classical and 
quantum capacity is calculated. The channel class appears in a study of the two-mode squeezing 
operator constructed from operators satisfying the fermionic algebra. We compare Grassmann 
channels with the channels induced by the bosonic two-mode squeezing operator. Among other 
results, we challenge the relevance of calculating entanglement measures to assess or compare the 
ability of bosonic and fermionic states to send quantum information to uniformly accelerated frames. 



I. INTRODUCTION 

The notion of quantum channel capacity is central to quantum Shannon theory. Early devel- 
opment in the seventies [l[ was a starting point to an impressive amount of knowledge that has 
been acquired in the last two decades @, Q ■ Two of the most investigated areas is the classi- 
cal [ij and quantum capacity Q of a quantum channel. The classical/quantum capacity informs 
us about the ability of a quantum channel to transmit classical or quantum correlations. More 
precisely, consider a sender who has, in principle, at his disposal an optimal encoder producing 
a classical or quantum code, and a receiver able to process the channel output and recover the 
transmitted information (that is, to decode) with arbitrarily high precision. In this way, the 
information can be transmitted at the rate given by the capacity and cannot be improved by 
any other choice of encoding. Various additional conditions or restrictions might be added, 
for instance if privacy is required [f| or some sort of assistance in terms of other quantum or 
classical resources available to the communicating parties 0, HI • This leads to a large number 
of important capacity definitions relevant under given circumstances and one might even try 
to characterize the whole capacity multi-dimensional regions in which the axes correspond to 
various available resources |9rfll|. 

Due to the presence of regularization (see below) the classical or quantum capacity is not effi- 
ciently computable. There are, however, particular examples of channels for which the classical 
or quantum capacity is easy to calculate. In the case of the classical capacity, every such exam- 
ple must be cherished since the proof usually involves some nontrivial manipulations fl2i - F]~7| . 
For the quantum capacity, almost all known non-trivial examples fall in the class of degradable 
channels [18|, |25|, |27[ . Among these examples are exceptional cases for which both capacities are 
known, to our knowledge there are only two examples: (i) a qubit erasure channel [l8| and (ii) 
Hadamard channels [14j |. There are also trivial examples of such channels with zero quantum 
capacity: entanglement-breaking channels [l5| and anti-degradable channels. In this paper we 
add another member into the elite group of non-trivial examples: the Grassmann channels. In 
order to find the quantum capacity we show that the Grassmann channels are degradable while 
to find the classical capacity we make use of the fact that the Grassmann channels are of a direct 
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sum form. We show that the channels which form the structure of the Grassmann channels are 
new members of the surprisingly broad family of channels studied in [ljj for which the classical 
capacity is efficiently calculable. The lowest-dimensional example of the Grassmann channels 
turns out to be a qubit erasure channel. In general, higher-dimensional Grassmann channels 
have certain traits in common with a qubit erasure channel making them interesting in the light 
of some recent capacity results [!, [l^| . 

This paper can be partially seen as an accompanying paper to Ref. [22]. There, one of the 
authors investigates a channel and its capacities induced by the bosonic squeezing transfor- 
mation, playing the role of a channel isometry. The isometry is presented from the physical 
point of view as it appears in the context of the Unruh effect for massless scalar fields. The 
channel was nicknamed the Unruh channel and its properties and analysis from the quantum 
Shannon theory point of view were also presented in other works EH [H, H3 ■ The cur- 
rent paper is the 'fermionic' version of the analysis done in the bosonic case, in the sense that 
the isometry is generated by operators obeying the canonical anticommutation relations (the 
fermionic algebra). The consequences for quantum Shannon theory, mentioned in the previous 
paragraph, are radically different from the bosonic case which justifies an in-depth study of this 
fermionic case. For example, we show that the Grassmann channels do not belong to the class 
of Hadamard channels, as opposed to the bosonic case [loj . In terms of physical consequences, 
the capacity calculations for the Grassmann channels are even more surprising. In contrary 
to common opinion, we show that: (i) it makes absolutely no qualitative difference for the 
study of the Unruh effect whether the states are composed of bosons or fermions (at least for 
quantum information transmission purposes) and (ii) naive calculations of certain entanglement 
measures of states shared between inertial and uniformly accelerated observers do not provide 
much insight into the channel's ability to reliably send quantum information. Finally, there 
is an intriguing connection of the Grassmann channels and their bosonic relatives to the well 
studied family of transpose-depolarizing channels. 

The encounter with fermionic degrees of freedom brings an interesting complication. Due 
to the use of fermionic statistics, the usual objects studied in quantum information theory, 
such as entangled or separable states, need to be treated carefully. As we discuss later in more 
detail, the problem lies in the fact that many-particle fermionic systems lack the tensor product 
structure. This is one of the simplest examples of the braided statistics and has relatively 
recently begun to be studied more closely by a number of authors [28| . Due to the presence 
of tensor products in all capacity definitions, it is not immediately clear how to generalize this 
concept to a fermionic system, or even a system with more general statistics. We will show that 
in our case, however, we may use the 'standard' framework for qubits (qudits) due to a careful 
choice of the input encoding for fermionic states and the specific isometry that we investigate in 
this work. Throughout this paper we use the multi-rail encoding and we will show that under 
given circumstances the capacity formulas indeed need not be modified. Thus, we may proceed 
and practice quantum Shannon theory with no modifications. 

The connection of quantum capacities and Grassmann variables has previously been briefly 
visited also in a different context by the authors of [29| elaborating on the formalism introduced 
in [H. 

The paper is organized as follows. In Sec. [TT] we recall some basic notions from quantum 
Shannon theory together with the fermionic algebra and briefly discuss various relevant physical 
situations. Sec. IIIII describes the construction of the qudit Grassmann channels and study some 
of their properties. Sections |TV] and [V] contain the calculation of the quantum and classical 
capacity of the whole (infinite-dimensional) class of Grassmann channels. In Sec. IVII we discuss 
several physical properties of the Grassmann channels. 
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II. DEFINITIONS AND NOTATION 



Classical and quantum capacity 

Let B(H) be the algebra of bounded operators on a Hilbert space H which for our purposes 
will be a d-dimensional complex vector space W — C d . Our algebra is therefore just the algebra 
of complex d-dimensional matrices. Let K, : B(Ha') ^ B(Ha') be a quantum channel mapping 
density operators from an input Hilbert space Ha' to an output Hilbert space T-La (if there is 
no chance of confusion we will use the shorthand JC : A 1 H> A). We define V/c : A 1 H> AC to 
be the channel's isometric extension. Equivalently, by embedding the isometry into a higher- 
dimensional Hilbert space we write U/c ■ A'C n- AC, where U/c is a unitary operator. Using 
the isometry picture the complementary channel of K, is defined as /C c (er) — Tr^fV/ccrV^] = 
Tta o V/c ° cr. The second equality is how we will occasionally abbreviate similar expressions. 
When there is no chance of confusion we will omit the mode index A or C for states. The 
von Neumann entropy H{qa) — — Tt:[qa log qa] will be written in the economic way as H(A) g 
which is mainly suitable for dealing with the parts of multipartite states qa = Ttbc... [qabc..]- 
This convention will also be used for other entropic quantities. 

A classical ensemble can be written in the form of a classical-quantum state axA' = 
J2xP( x )\ x )( x \x ® Vx,A' with a probability distribution function p(x). Then, for a quantum 
channel K, the capacity is given by the regularized expression Q 

C{K) = lim -C^{K® n ). (1) 

n— too XI 

Choi(^C) is the Holevo quantity 

C Ho i(/C)= sup I(X:A) q , (2) 

{P( X ),<?X,A>} 

where <,xa = fc{&XA') and I(X : A) = H(X) + H(A) — H{XA) is the mutual information. For 
the optimization task in Eq. ([2]) it is sufficient to consider <j x .a' to be pure states. 
The expression for the quantum capacity contains regularization as well [B[ 

Q(/C) = lim -Q«(/C®"). (3) 

71— S-OO TL 

Q^(JC) is the optimized coherent information 

QV>(lC)=sap[H(A) T -H(C) T ], (4) 

</> 

where t ac = Vk, ° ipA> ■ 

There exists a generalization of the von Neumann entropy known as the a-Renyi entropy 

H a (p) = 1/(1 - a)log[Trp Q ], where the von Neumann entropy is recovered for a — > 1. A 
quantum channel V is called degradable [25|, [H[ if there exists another channel M. such that 
A4 o T) = T> c holds. The channel A4 is called a degrading map. Finally In is the natural 
logarithm and log denotes the base d logarithm unless stated otherwise. The motivation for 
using the base d logarithm is a better way of comparing the capacities for the class of channels 
studied here and the their 'bosonic' equivalent studied elsewhere (2^ . 



Fermions and their correlated pairs 



For a set of modes that are specified by quantum numbers compactly denoted by k, the 
exchange characteristics of indistinguishable fermions that may occupy those modes are reflected 
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in the anticommutation relations obeyed by the corresponding creation (annihilation) operators 

a « (4) : 

{a K ,al,} = 6{k- k'), (5) 
{«L a l'} = Wn,a K '} = 0. (6) 

An interesting transformation that preserves the anticommutation relations results from the fol- 
lowing algorithm. Within the available modes two subsets are chosen. Each of the modes from 
the subset of lower cardinality is related through an injective function /(k) = k' to the elements 
of the other subset. In this way, the modes are paired and so are the creation-annihilation op- 
erators a K , 4 and ct/(«), a /( K )- I R order to make simpler the notation, the operators a/( K ), a /( K ) 
from now on will be denoted by a different letter, for instance, a j ^ , a J, / fi s — > c K ,cJ.. It can be 
directly shown that the so called Bogoliubov transformations 

b K — cos ra K — e _i * sin rc\ (7a) 

b\ = cos raj. — sin rc K (7b) 

d K = cos rc K + e -1 ^ sin raj. (7c) 

<4 = cos rcj. + e 4 * sin ra K . (7d) 

preserve the anticommutation relations Eq. ([5]). These transformations were first introduced 
for the generation of the Bardeen-Cooper-Schreiffer (BCS) states, which are an excellent ap- 
proximation to the ground state of a weakly interacting superconductor [32| . In that case, the 
modes are thought to describe electrons with the quantum numbers usually taken as the vector 
wave number k and the projection of the spin s. The function / is such that f(k, s) = {—k, —s). 
The Cooper pairs that are described using the operators b R , d K , as a consequence, have opposite 
momenta and spin. 

Bogoliubov transformations have been very useful in the description of strongly correlated 
fermions in diverse scenarios like condensed matter [33j |. ultracold degenerate atomic Fermi 
gases [HI or quantum field theory of particles on spacetime with a nontrivial metric [H, H(| 
|H, H^l . To illustrate the impact of these transformations in the latter area, consider negligible 
interacting massive fermions. For a uniformly accelerated observer, the adequate spacetime 
coordinates are the Rindler ones; for which two wedges, that can not be causally connected, are 
identified, we shall denote them by right R and left L wedges. The right (left) Rindler modes are 
only supported on the right (left) wedge. An important question concerns to the energy spectra 
seen by a Rindler observer in connection to the Unruh effect [13] ■ The most clear treatments of 
the problem rely on comparisons between the so-called Unruh modes and the Rindler modes. 
They are naturally carried out in terms of generalized Bogoliubov transformations [35| that 
take into account that each Unruh mode needs for its representation an infinite superposition 
of Rindler modes. In this kind of scenario, perhaps the simplest realization of a Bogoliubov 
transformation corresponds to the case of a spin one-half fcrmion moving in a one dimensional 
space. Then, the modes of each fermion can be described in terms of a Grassmann field and a 
basic transformation is 

( «M \ _( cosr -e-^sinA ( b k , s \ 

Ulf-J" V^sinr cosr ) \P k ,_ J ' (8) 

where a, stand for particles and a, a) for antiparticles, k denotes the wave number and s its 
spin projection. The real parameter r can be chosen to depend on the acceleration a and the 

rest mass m of the Rindler observer as tanr = e^ 2 ^ \/ k2 +( mc / h ) 2 [38] . Note that Eq. ([SJ is 
Eqs. (17a[) and (|7dp written in the matrix form and with the notation adjusted to the relevant 
physical situation. 
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III. GRASSMANN CHANNELS 



Definition 1. We define a d-mode fermionic state as 

d 



\ F ) = Il^lr |«oc) = \nx . . . n d ) = \n), 



where rii € {0, 1}. Following the properties of the fermionic operators we see that (i) each mode 
is occupied by at most one particle and (ii) the state Ji* 1 ) is completely antisymmetric. For 

< k < d there is (^) possible fermionic states for which J3 i=1 ni = k holds. 

Note that throughout this article all input Hilbert spaces are spanned by the multi-rail basis. 

Definition 2. The multi-rail basis is defined as the set of all states \F} for which ni = 1. 

The basis spans a d- dimensional Hilbert space and any fermionic state written in this basis is 
said to utilize the multi-rail encoding. 

The operator related to a Bogoliubov transformation Eq. (|8]) over the i-th fermionic modes 
a and c reads 

U AiCi = exp [r(a\4e-** - c t a^)} , (9) 
where r,cj> € R- The operator exponent may be factorized according to the following theorem. 
Theorem 3 ([21]). Let J + ,J_ and J3 be operators satisfying the commutation relations 

[J 3 ,J±] = ±J± 
[,/+,J_] = 2J 3 . 

Then the following identity holds 

e X + .7 + + X 3 .7 3 + X-J- _ e A+J+ e lnA 3 ./3 e A- J- ^ £q\ 



(10) 



where 



2A± sinh f 

A+ = 



2/ cosh / — A3 sinh / 
A 3 = (cosh/- ^ sinh/) ( 12 ) 
/ = ((A 3 /2) 2 + A_A + ) 1/2 . 



In our case we have 



J+ = aV 

J- = ca = —ac 

■h = \ (a f a + c f c - 1) 



(13) 



This choice satisfies the commutation relations in Eqs. (fT0|) . We get from Eq. Q A± = ire^ 4 * 
and A3 = leading to 

UAid — cosr exp \a\c\e~ % ^ tanr] exp [(ajaj + cjcj) In cos^ 1 r] exp [— Cjaje 1 ^ tan r]. (14) 

We might safely set = since we will later see that it has no relevance in this work. 

We collect the majority of identities used in the course of the paper in the following lemma. 
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Lemma 4. Let a and c be operators obeying the canonical anticommutation relations Eq. f5|). 
Then the following identities hold 

(15a) 
(15b) 
(15c) 
(15d) 



(15e) 
(15f) 
(15g) 



The first equation is a generalization of the second row in Eq. UU\), Eq. \15c\l simplifies to 
t„.„t _ „t 



a i a i a i = a i f or ' = 3> Eq. H5e]) holds for i ^ j and Ak-i — k(k — l)/2. 

Proof. Eqs. (|15a|) - (|15fj) directly follow from Eq. ©. In Eq. Q15gD we move all dj on the left of 
all Cj. So all dj on the LHS of Eq. ( |15g[ ) for even j 'jump over' an odd number of Cj operators. 
Similarly every odd aj switches its position with an even number of Cj operators. The total 
acquired phase is 



[oiCi,atct 


1 = 25 4 A 


[a\a u a^c] 


] = ^-atcj 




1 = (S^-a 1 
J u y "j 


[a\o),a\ 


] = o 




- k 








= 




-3=1 














= 




i j=i 




k 


k k 


(") 


Afci n^ c 




3=1 





TT(_)(fe-j) = (_)E*=i(fc-i) = (_)*(fc-i)/2 = (-) A *<-i. 



Taking into account Eq. Q15g| the action of d copies of the fermionic unitary operator results 
in 



I*) 



AC 



\Ua z c, \vac) = cos d r y^(-) Afc - 1 tan fc r ^ \™)a 



(16) 



fc=0 



If we wanted to see how Uac = ®i=i U Aid transforms a fermionic qudit written in the 
multi-rail basis 



\^)a'c = ^2 \ vac ) = X! & V)a> \ vac ) C " 



we might just calculate = Uac I^a'c ^ n ora ^ er to simplify this complicated calculation 

we first observe 

d d d 

Uac — cos d r exp tan r^^ajcj exp |^y^(a|ai + c\cj) In cos -1 r exp tan r^^a^ , (17) 



i=l 



»=i 
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where Eqs. (|15a|) and (Il5b[) for i ^ j were utilized. To proceed we make use of the main 
advantage of the multi-rail encoding. Due to the presence of Ci annihilating the vacuum state 
we observe 



exp 



a a a 

tanr^^aiCiJ ( A a lj \vac) = /^aj \vac) 



i=l 



We may then write 



a a a 

\®) AC — cos d r exp tanr^ajct exp |^E (aja; + etc*) In cos -1 rj (E ftaj) |uac) 



cos d 1 r exp 



tanr^ajctj (E&°<) \vac) 

d d 

cos^^f Vftajjexp tanrE a ! c ! |wac) . 



(18) 



The second row follows from Eq. (|15cl) (see why i = j does not spoil the commutator). The 
last equality is possible due to Eq. (|15d[) . 
The action of the unitary Uab leads to 



\$) Ac = U~Ac'^2fii\i) A , \vac) c , = cos d 1 r ( E ftaj j exp tan^E a]c\ \vac) 



3=1 



d+1 



»=1 fc=l JV fc 



d 



= cos d 1 r E tan?i: 1 r (- 



fc=i 



d-fc+l 



E E (±)«a 



L JVfc iei 



(19a) 
(19b) 



(19c) 



where we have defined the set Nk = {(fii, . . . , rid) I X) 7 =i ' 



1} to be the sum over all 



possible states with k — 1 fcrmions. The additional fermion (with respect to the C subsystem) 
in Eq. (|19c[) occupies the i-th position. There is states | . . . rij . . . ) c in the middle sum 

of Eq. (I19c|) for a given jVfc (fc = 1 . . . d). The rightmost sum of Eq. (|19c|) sums over i € I, for a 
fixed JVfc, such that / = {i | ni = 0} and thus contains rf— fc + 1 terms, justifying the upper limit 
in the sum. Taking, for example, flo = 1 we verify cos 2 ^ 1 ) r X)fe=i tan 2 ^ -1 ) r(tZi) = 1. The 
multiplicative sign (±)j in Eq. (|19cl) is dependent on the fermionic state of the A subsystem. It 
appears following the operator rules specific to the fermionic algebra. Let us recall some basic 
properties of the algebra 



a\\...rij...) = (l-n;)(-) 
Oi| ...rij...) =n l (-)^=i" 3 



12) Jl n 3 



. Tlj + l. t 



> 



(20) 



Isometry output Eq. (|19[) gives rise to a new class of channels we nicknamed Grassmann chan- 
nels. 
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Definition 5. Let the d- dimensional Grassmann channel Gd '■ B(C d ) >— > B(C 2 _1 ) be a quantum 
channel defined by the action of its isometry Vg d as Gd(i>A') = Trc o Vg d o i/ja' ■ The action of 
the channel on an input audit is given by 



'd-r 

,k - 1, 

fe=l v ' k=l 



G d ■■ 4>A' H« if = cos 2 ^ 1 ' r 0tan 2 ( fc - 1 ) _ \) X [ d) = © Pfc &,*(iM> (21) 



w/iere p fc = cos 2 ^ 1 ) r tan 2 ^ 1 ) r(fz\) ■ The maps Gd,k ■ B{C d ) ^ fi(c(*)) for 1 < k < d 
are quantum channels constituting the Grassmann channel Gd and Xk = Sdki^A') having 
Tr[xi d) ]=L 

Remark. The output block dimension (f\ coincides with the dimension of the vector spaces into 
which the Grassmann algebra A(V) over a d-dimensional vector space V decomposes (A(V) = 

0fe=l Afe ( F ))> hence the 

name Grassmann channels. Note that the Grassmann algebra is also 
known as the exterior algebra. 

Remark. The index k in Eq. (|2"Tj) has two roles in \ f ■ It labels the state but it also indicates 

how many fermions the basis is composed of. We thus know in which basis xjf^ * s written. 

Remark. The first Grassmann channel Gi is a trivial trace map. The second Grassmann channel 
G2 is a channel recently playing an important role in quantum Shannon theory - a qubit erasure 
channel [3j, Il8l420l | . We now look at Gi and G3 cases in more detail. 

Example. For d = 2 we get from Eq. (Qi 



cos 2 r 



cos 2 ' 



|^) AC . = cos 2 r exp tan r(a[c[ + a 2 c 2 ) \vac) 

1 t t t t\ tan 2 r. 4 + tt\2 1 \ 
- tan r(a[c[ + a 2 c 2 ) H — {a[c[ + a' 2 c 2 ) \vac) ( 2 2) 

' r(^\vac)+ tan r( |l0>^ |l0> c +|0l)^ |0l) o ) - tan 2 r |ll) A |ll) c V 

Following Eq. (|19cp we obtain 

\$) AC = cos r ((ft |10) +/3 2 |01» |uac) + tanr (ft |ll) |0l) -ft |ll> 1 10» ) (23a) 
= ^/r^(ft|l>+ft|2>) A |0> c +Vp|3> A (ft|2>-ft|l>) c . (23b) 



r 



The last equation is Eq. (|23al) rewritten using a reparametrization -^/l — p — cos r, ^Jp = sin 
and a logical ket notation (|0), . . . , |3)) to facilitate the comparison with an isometry output 
for a qubit erasure channel 



WIT = V^(/3i |1> +ft |2> ) A \f) c +VP\I) A (A |1> +/^2 |2) ) c . (24) 

J /) is a flag state orthogonal to both 1 1) and 1 2) . The most notable difference between Eq. (I23bl) 
and Eq. (|24]l is the dimension of the output Hilbert space. For the latter the Hilbert space 
is three-dimensional. But if we trace over the A or C subsystems and compare them we 
immediately see that they indeed induce the same channel. The Grassmann channel G2 is just 
embedded in a higher-dimensional space than the corresponding erasure channel and unitarily 
rotated (note the difference in the second brackets of Eqs. (|23bj) and (j24jl ). 
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Example. For Q3 we get from Eq. (|19cl) 

|$) Ac = cos 2 r((ft 1 100) +ft |010) +ft |001» |«ac) (25a) 

+ tanr [(-ft |l0l) -ft |011» |00l) -(ft |ll0) -ft |011» |010) +(ft |H0) +ft |101» |l00>] 

(25b) 

-tan 2 r |lll)(ft |01l) -ft |l0l) +ft |110>)Y (25c) 

The trace over the C subsystem gives us the output of G3 

£3 : $A> ^ if = Pixf © P2X ( 2 ] © 2>3x! 3) , (26) 

where pi = cos 4 r,p 2 = 2 tan 2 r cos 4 r, p 3 = tan 4 r cos 4 r. xf is the input state itself and xf 
is a flag state We are interested in the form of x^ 

1 /lft| 2 + lft| 2 ftft -ftft 
Xa } = J &ft lft| 2 _+lft| 2 ftft I (27) 



2 



-ftft ftft |/3i| 2 + |ft| 2 



where the two-fermionic basis is ordered as {|01l), |l0l), |ll0)}. We easily verify that the 
Kraus operators 

, /o -1 o\ 1 /o i\ , /o o> 

#! = -= ! ,if 2 = — \,K 3 = — 1 I (28) 



^ \0 0/ ^100/ 



v0 1 0; 



yields xf' '■ 

Remark. The basis in which the states and Kraus operators from the two previous examples 
are written are generated by the action of creation operators satisfying the canonical anticom- 
mutation relations. Yet in our case we treat the Grassmann channel outputs as if written in 
a multi-qubit basis and we apply the usual machinery of quantum Shannon theory (namely 
utilizing the tensor product structure). The following Theorem justifies this step. Note that we 
will illustrate it on the case of the quantum capacity. The classical capacity follows the same 
line of arguments. 

Theorem 6. The quantum capacity formula Eq. |3J) remains unchanged if for an isometry 
Vk. ■ A' AC the input Hilbert space is spanned by the fermionic multi-rail basis introduced in 
Definition [H 

Proof. Using the definition of the optimized coherent information Eq. Q we write the relevant 
part of Eq. d3j) for n = 2 



(/C® 2 ) = sup [H(AxA 2 ) T ^ - H{dC 2 ) T ml (29) 
e 



where O = Qa^C A'^C 2 is an input state. The structure of the input state is crucial. First, 
recall from Sec. [n] that we represent isometry Vjc as a unitary from Eq. ([17} U ac '■ A'C H> 
AC where the input reference system C is prepared in a vacuum state. Therefore, every 
pair of input modes A'jC'j contains exactly one fermion and an arbitrary input basis state 

{ 1} A'C A'C }fe— 1 con t ams two fermions. The basis is clearly constructed by tensoring 
product basis of two Hilbert spaces A' 1 C[ and A' 2 C' 2 . Note the order of the product reflected by 
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the order of the primed subscripts. We now turn our attention to the action of (U A( L ® ^ac) 
on this basis. The superscripts on the unitaries go from one to n to distinguish it from J/a^c* 
in Eq. (JTB]). We find that 

® Uf c ) \^) A , C , A , C , = \^) AC ® |*^> AC - r« CiA2C2 , (30) 

where the middle equation is a product of two states from Eq. (|19cp each of them having all /3's 
but one equal zero (fci, k 2 = 1 . . . d). This crucial step is possible due to the fact that each pair 
of modes A'jC'j contains just one fermion, the structure of U AC in Eq. (jTTj) and identity (| 15d[) . 

Looking at Eq. (j!9cj) (or better at the examples in Eqs. (I23a[) and ([231) representing a generic 
situation) we see that 

SWAP AaCl r^ ClAaCa * r A k) iMCiC2 . (31) 

A sign change occurs according to the number of fermions the swapped systems contain. How- 
ever, when we apply the unitary product on an arbitrary input state 

{uTc®uf c )l^a k \Q^) A , c , A , c ^ =r AlClA2C2 , (32) 

where |afc| 2 = 1, and reorder the output followed by partial trace over C\C 2 (AiA?) to get 
the (complementary) channel output, we notice another crucial property. The partial trace 
produces a direct sum structure where for each subspace the sign that comes from the reordering 
procedure is identical. The reason is that each block contains only fermionic states of a given 
fermionic number and its constituents come from the parts of t Ai c 1 a 2 c 2 where the parity of 
the number of fermions in C\ and A2 remains unchanged (so the 'swapping' sign is the same). 
So if the sign is positive (even or odd number of fermions in both subsystems) nothing happens 
and if the sign is negative (opposite parities in C\ and A 2 ) the sign gets canceled for a given 
subspace. Stated differently, the phase (in this case just plus or minus one) is common for each 
subspace so it disappears in the density matrix formalism (blockwise) and therefore for the 
whole block-diagonal density matrix. 

Only little changes if we take a tensor product of n unitaries in Eq. (|29[) . The dimension of 
the input Hilbert space will be d n and the equivalent of Eq. (|30p will hold as well as the phase 
argument for the resulting block-diagonal structure after taking partial trace over the C\ . . . C n 
or A\ . . . A n subsystem. ■ 

Remark. Consequently, provided that we use the multi-rail encoding for the isometry input, 
we may simply take an output state Ta(C) an d pretend that the basis in which it is written is 
an ordinary multi-qubit basis. Thus taking the tensor product in quantities defined in Eqs. ([1]) 
and ([3]) is justified by the above theorem. We may take their tensor products following the 
usual rules as in ordinary quantum information theory. 

Also note that in the proof of degradability of Grassmann channels in Section IIVI we work 
with entire isometry outputs and so we have to fully respect the fermionic character of these 
states. 

In the proof of the next theorem we will extensively use the geometric picture from the rep- 
resentation theory of the s\(d, C) Lie algebras. We summarized some basic facts in Appendix [A] 

Theorem 7. Let the first block of 7^ in Eq. \21)) be written as 

xf = ~(t + J2n a \U), (33) 

ct=l 
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where \^ are generators of the dual representation to the fundamental representation of the 
s\(d, C) algebra, L = 2d(d — 1) and n a are functions of fiiftj- Then the remaining blocks in 
Eq. \21}) can be expanded with the same coefficients n a 



(d) 



a=l 



(34) 



where are generators of the k-th completely antisymmetric representation of the s\(d,C) 
algebra. 

Remark. The number L in Eqs. (|33f) and Q34p comes from the use of a redundant number of the 
sl(2, C) subalgebras as discussed in Appendix [5] Namely, there is d(d — 1) diagonal generators 
H n , d(d— l)/2 off-diagonal generators E^j (recall that i < j) and d(d— l)/2 of their Hermitian 
conjugates e\-. 

Remark. Since we use a linearly dependent set of algebra generators the expansion coefficients 
n a are not unambiguously determined. For us, however, it is sufficient to show that at least 
for one specific construction (the one presented here) the coefficients can be chosen to stay 
preserved when switching to a higher-dimensional representation of si (of, C). 

The upcoming proof will be divided into two parts. We will separately prove Eq. (|34[) for 
diagonal and off-diagonal generators. To better follow the proof it might be helpful to watch 
the example of Eq. p9c[) for d = 4. In the course of the proof there is a remark illustrating 
results on this case. 



Example. For d — 4 we get from Eq. (|19c 



\$) AC = cos 3 r^^i |1000)+Aj |0100)+/3 3 |0010)+/3 4 |000l)) \vac) 
+ tanr (/3i |l001>+/3 2 |0101>+/3 3 |001l)) |000l) 



1001)- 


f/?2 


0101) 


+/3 3 


0011} ) 


0001) 


1010)- 




0110) 


-/3 4 


0011)") 


0010) 


1100) 


-ft 


0110) 


-ft 


0101)") 


0100) 


1100)- 




1010) 


fft 


1001)") 


1000) 



tan 2 r 


(jSi|1011>- 


hft |0111>) 


0011)4 


(ft 110l)-ft 0111)) 


0101) 




(ft 1101)- 


fft|l011>) 


1001)4 


-(ft 1110) +/3 4 0111)) 


0110) 



(35) 



- tan 3 ; 



1 1111) (-A |011l)+ft |l01l)-ft |ll0l)+ft |1110) N 



!>) I 1100 ) 



We do not indicate the subsystems but in a product of two kets the first one is the A subsystem 
and the second one is the C subsystem. 
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We will employ the fermionic representation of the si (2, C) algebra (see also Eqs. (|A3[i 1 

(36a) 
(36b) 



1 

-1 



= a] n a m - ajai 



1 





= oLai 





1 



a a 



(36c) 



where m ^ I are mode labels. 



Proof, (i) Let us first find the coefficients n a from Eq. (|33p and then show that they remain 
preserved in Eq. (f34|). Since we have 



1/3- 12 ' 



d-l 



iL( 5>}o, + - a t 1+J a 1+J )) = |ft| 2 aki = |/3 X | 2 |1 . . . 0)(1 . . . 0| 



(37) 



and so 



d-l 



^f- ( ^ ajai + J2( a U* ~ «lei^©i)) = IftlV* = I AI 2 I° ■ ■ • ■ ■ • °X° • • • • • • °l ( 38 ) 

i=l j=l 

there is (d— 1) coefficients n Q that equal |/3i| 2 for i = 1 . . . d. Note that i®j = i+j mod (d+1). 

For fc > 1 we first notice in Eq. (I19c[) the presence of the factors (±)j. They are irrelevant 
for now since we investigate the diagonal generators. The important fact is that the numerical 
coefficients of all fa's are the same in the absolute value, that is, equal to one. We trace over 
C and the diagonal part of the A subsystem for a given k reads 



diag[xi d) ] = -^y 



1 



d-k+l 

]T ]T |A| 2 |...n J - + l i 
L N k iei 



. Uj + 1 



Each diagonal element Yltei +l \Pi\ 2 \ ■ ■ ■ n j + 1» 
can be factorized 



(39) 

corresponding to a given rij 



d-k+l 

\(3 l \ 2 \...n J +l l ...)(...n J + l l ...\ 

iei 

(d-k+l 
Ift| 2 |0...1,:...0)(0...1 t ...0| 



"3 • ■ • I ) 



(40) 



where the expression in the parenthesis in the second line can be rewritten using the fermionic 
representation of the si (2, C) algebra from Eq. (|38|) . The rest serves as a label (recall that there 
is (t_i) orthogonal states for which ^2 n j = k — 1). However, because of the factorization in 
Eq. (|40p the dimension of the first expression will become d' = d — 5~]j—i nj giving us 

X- • • nj + U . . . | 



d-k+l 

£ |A| 2 |...n,- + l i ... 
iei 

/d-k+l d' d'—l 

= ( \Pi\ 2 (Y a \ a i + ^2( a l a i ~ a i®j a iBj) 



iei 



■)(■ 



(41) 
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So the diagonal elements for k > 1 can also be expressed using only the sl(2, C) diagonal algebra 
generators. This is what we expected following the discussion in Appendix [3] 

Now we have to make sure that there is the right number of summands when constructing 
the direct sum sl(2,C) subalgebra representations also discussed in Appendix lAl in the last 
paragraph. This is indeed the case. The orthogonal 'label states' factorized out in Eq. (|4"Tj) 
label the subspace in which the si (2, C) subalgebra in the parenthesis lives. We want to count 
how often this situation happens. Since it is a si (2, C) subalgebra we have two spots out of d 
in each ket occupied by two fermions (by choosing, for example, /3i and fa or any other pair) 
and the remaining number of spots can be occupied by n j fermions. The number of 

possibilities is 

^[orthogonal subspaces] = ( 

VEj 1 11 j 

We thus reproduced the result found in Appendix [A] for all completely antisymmetric represen- 
tations of the s\(d, C) algebras. 

Remark. We might verify the above considerations on Eq. (|3"5|) . For k = 2 (lines 2 to 5) tracing 
over C and taking the diagonal part gives us unnormalized Eq. (l39l) for which d — k + 1 = 3. 
There is (^) = 4 states where J2j=i n j = 1 holds and they can be factorized out as in Eq. (|4"Tj) . 
There are therefore four blocks with mutually orthogonal flags attached to them but only (^) = 2 
of them have a common pair of coefficients (e.g., fa and fa). This is the case illustrated in 
Fig. [SJd. All edges of the octahedron are sl(2, C) algebras. The pairs of parallel lines form the 
generators of a completely antisymmetric representation of si (4, C) in the form of a direct sum. 
They span a 4-dimensional space but none of the spaces is orthogonal to any other. 

(ii) Proving Eq. (|3~4"|) for off-diagonal generators is considerably simpler. Looking at 
the rightmost sum of Eq. (|19c|) we see that every state | . . . rij . . . ) c is multiplied by 
Y]j f3j(±)j j . . . rij + lj . . . )^ containing one fermion more in the i-th mode compared to the 
C subsystem. Tracing over C and taking the off-diagonal part gives us expressions of the 
following type 

A£m(±M±)m| • ■ • m + h ...)(... nj + l m ... |. (42) 

l.rn 

We immediately recognize the off-diagonal matrices forming Eq. (|42[) to be proportional to 
the step operators pertaining to the sl(2,C) subalgebra from Eq (|36|) . The situation is only 
complicated by the presence of a function (±)z(±) m responsible for the change of a sign. The 
function comes from the fermionic relations Eq. (|2T)]) . The change of sign therefore appears if 
we annihilate the fermion in the l-th mode and create it in the m-th mode 

(±)i(±) m j . . . n 3 + l m . . . ) A = a^ai \ . . .rij + li . . . ) . 

But this is precisely the action of the operator from Eq. ([36bJ representing the sl(2, C) subalge- 
bra. We again take a direct sum of ( k Z^) off-diagonal generators. They form the off-diagonal 
generators of the A:-th lowest completely antisymmetric representation of si (d, C) . H 

Corollary 8. The Grassmann channels Gd are SU(d) covariant. 



Proof. The proof is identical to Corollary 17 for completely symmetric representations of sl(c?, C) 
in [22] showing that the Unruh channel is SU(d) covariant. ■ 
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IV. QUANTUM CAPACITY 

There is an interesting symmetry between the A and C output subsystem captured in the 
following lemmas and later necessary for the proof of degradability of the Grassmann channels. 

Lemma 9. Labeling the Pauli matrices o~z cind ax by Z and X and a d- dimensional identity 
by Id we first introduce an infinite product Z®l2®Z®l2®Z®l2®.... A d-mode unitary 
operator is defined as the first d unitaries of the product. Then the following identity holds 

d d 

Oo-Kd-i) ° X ® 2d ® Id) ° H d (Eft4) \vac) AC = tan^ 1 r( ]T ftc]) \n) A \n) c , (43) 

<=1 i=l 

where X)i=i n * = ^ — 1 an< ^ C(fe-i)->(d-fc) : tan fc_1 r — > tan d ~ fe r is a non-physical operation solely 
acting on the trigonometric function in Eq. \19b\) (or Eq. H19c\) ). For k = 1 we understand 

O -Kd-i)(l) = tand_1 r m Ea - <Ef- 

Proof. We prove the identity by a direct calculation. Using properties of the fermionic algebra, 
the following hold: 

d d 

X® 2d o (E d ®I d )(22 ftaj |«ac) AO ) = X® 2d o (S d ® 7 rf ) o ft |0 . . . U . . . 0) A ) |«ac) c 
t=i t=i 

(44a) 

d 

= { -y X ® 2d o ft |0 . . . U . . . 0) A ) |mc> c (44b) 

z=l 

d 

= (-) i (EA|l...O i ...l) A ) |1...1> C . (44c) 



On the other hand we have 
d 

tan 

i=l 



a a 

d -' r A4 \n) A \n) c = tan"" 1 r ( ^R^+^ft |l . . . 0* . . . l) A ) |l . . . 1> C . (45) 



By multiplying Eq. (|44al) by (— ) d 2 and exchanging the trigonometric function (represented by 
the action of Oo-^(d-i)) we g e ^ the RHS of Eq. (|45|) and therefore also Eq. ((43]). ■ 



The purpose of the next lemma is to present two extensively used identities. 
Lemma 10. Let a be a fermionic operator. Then, the following identities hold 

Z<J = -a)Z (46a) 
Xa f = aX. (46b) 

Proof. Directly follows from Eq. §5§ once we realize that X = a + a' and Z = a^a — aa) — 
2at a - 1 = 1 - 2aa t . ■ 

Lemma 11. The following relation holds 

d d 

CVu-Kd-*) o X® 2d o (E d ®I d )o(jT ftaj) \*) AC = ( 04) I (±U,k,d*) AC , (47) 

i=l i=l 

where | (±) 7 ?i,fc,d^ r } AC is \^>) AC from Eq. \16\) with some of the summands having a negative 
sign. 
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Remark. The indices of the phase function will be described and the function explicitly written. 

Proof. To make the derivation smoother we will first prove a specific part of Eq. (1471) in which 

k = 2 



a 

Oi_>(d_ 2 )(taar) ° X md o (E d <g> / d ) o (^ftafjatej l mc ) 



AC 



i=l 



tan^^-y-^^Acj) I^Wc- (48) 



i=l 



where TO i = ^ ~ 2 and = 0. Let us ignore for a while the function Oi_>(d_2)- Using 

identities (I15dp . (I46al) . (l46b[) and Lemma [5] (note the missing j-th summand from both sides of 
Eq. @3)) the LHS of Eq. g8]) reads 



a a 



i=l i=l 



a 



(49b) 



where T = (— 2 (— ) J 2 = (— ) J 1 (— )■? 1 = 1 depending on whether i < j or j < i. Thus 
we recovered the phase on the RHS of Eq. (|48|) . It is time to justify the presence of Oi^(d-2) 
in Eq. (|4"5)l . Eq. (I19c[) indicates that the exponent of tan'""" 1 r coincides with the number of 
excitations of the C subsystem and is one less than the number of excitations in the A subsystem. 
From reasons that will become clear later we would like the exponent of the trigonometric 
function tan d_2 r on the RHS to correspond to the number of excitations of the A subsystem 
which is indeed equal to d — 2. 

To prove Eq. (|4T|) we first rewrite Eq. (IT5|) as 



d+l d+i 1 

£ tan- 1 r(-)*- £ \r?) A \H J ) C = £ tan*" 1 r 5>J c t)»i = __ j*)^, 



fe=i Wfc fc=i Wfc 

where we have used the notation introduced in Eq. (|19c[) . We have also introduced the set J, 
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for a fixed Nk, as J = < j < d, rij = 1}. Let us investigate the expression for a chosen k: 

d 

(fe -i H ( d - fe )(tan fe -V)oX® 2rf o(H d ®7 d )o(^^a])^;(at c ])^ |vac) AO (50a) 



i=l N k 



(fc -i)^ ((i ^)(tan fe - 1 r)o^((-l)^ J c j r j o X^ 2d o (S d ® 7 d ) o (^ftajj |«ac) AC (50b) 

d 

tan d - fe r(-)^((-iy ajCj )^ (£>ej) |rT') A |n') c (50c) 

tan d - fe r(-)' i ^(-)^(-)(' £ - 1 )( d - 1 )T' £ - 1 (-) A -^^/34) |m J ) A |m J ) c (50d) 

d 

tan d - fe r ( g ftct) ^(±) A , M (m^ |m J ) c , (50e) 



i=l JV fe 
i£7 



In the first equality of Eq. (1501) we used Eqs. (1 1 5 f f) . (|46al) and (|46bp . and in the second equality 
Lemma [9] was used with the summands corresponding to the set J removed from both sides of 
Eq. (|43j). Here J^i n' i = d—1 holds. In the third equality, Eqs. (I15c|) . ( |15g[ ) and the properties 
of the fermionic algebra leading to T in Eq. (|49b[) were used and we introduced a vector \rh J } 
such that X)j=i m f = d — k and to^ 7 = if j e J. Finally in the last equation we collected all 
phases into a single function (±),n,k,d- 

To prove the lemma we realize that the LHS of Eq. (|47|) is Eq. (|50ap summed over /c 



cos r 

fe=i 



•E 



0( fc - 1) ^ (d - fe) (tan^ 1 r) o X« 2d o (E d ® J d ) o (V ftoj) £(o}ct) 



i=l iV fe 



>ac ( 51 ) 



(note that we sum only to c? since the (d + l)-th summand disappeared). The RHS of Eq. (l4"Tl) 
is just a sum of Eq. (|50e|) over k 

d d d 

cos d r^tan d - fe r(^/34)^(±) A , M |m J ) A |m J ) c = |(±).a,M*>AC ( 52 ) 

fe=l i=l JV), i=l 

Notice that if we exchange (d — fc) for (k — 1) in ta.n d ~ k r, Nk, 4 an d \ T ™ J )c sucn that 
Ei=i TO / = fc — 1 we indeed get ( ^ ftcf ) l^)^ up to the phase function. ■ 



Remark. The important fact about Eq. (|50cl) is that the sign depends on chosen d, k and Nk 
so it is in general different for each summand over Nk but common for each sum over i. 

Lemma 12. Following the notation of Definition^ we label 7^ = Tr^ o <&ac- Then we have 

{k -i)^{d-k) ° W o 7 ^ = ^ , (53) 
where W is a unitary transformation. 
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Proof. Lemma [TT] is powerful since it allows us to compare the outputs of Grassmann channels 
with their complementary outputs 7^. Indeed, from the form of |3>)^ c we can deduce the 
explicit form of the block matrices from which 7^ is composed of (see Eq. (|2ip). However, if 
we wanted to compare 7^ with 7^ given by unitary U ac hi Eqs. (|19[) it would be a difficult 

task. The advantage of Lemma [TT] is that we don't even need to know 7^ or 7^ explicitly to 
find the relation between them. 

Let us rewrite the previous lemma result in the following way 



d d 



\*) AC = (k _ 1)Md _ k) o(~ d ®I d )oX® 2d o ($>cj) |(±U M *> AC . (54) 



i=i i=i 

Recall that 0( k _i- ) ^( d _ k - ) o(E. d ®I d )oX® 2d is an involution and 0(fc-i)->(<j-fc) is a scalar function. 
Note that S d <8> Id commutes or anticommutes with X® 2d depending on the specific form of S^. 
We incorporate this sign change into {^)rh,k,d- Tracing out the A subsystem we have 

d d 

7 £°=TV A o(£>aj) \V) AC = (k _ 1) ^ (d _ k) oW oTt a o(J2 Pi4) | (±)tS,*,<i*> ac (55a) 

Z=l 1=1 

d 



(fc _ w) o^oTV A ogft c } \^) Ac (55b) 



j=i 



= 0(fe- 1 )^(d- fe )°WoTr c o(^ftatJ |*> AC (55c) 

i=i 

= 0( fc -i)^( d - fe )°^°7i d) - (55d) 

The first equality holds since all the unitaries in Eq. (|54[) act locally and O (u_x\^i d _f.-s is again a 
scalar function (the unitary W appears as the result of tracing over unitarily- locally transformed 
state The sign ambivalence in the second equality is irrelevant because tracing over A 

means creating of a convex sum of states belonging to the C subsystem so the phase disappears, 
the third equality is a simple permutation argument due to the symmetry of |^') j 4 C . between 

the modes A and C and in Eq. (|55d[) we invoke the definition of j A . ■ 

Let us formally define the complementary Grassmann channel based on Lemma [T^] and Def- 
inition [5] 

Definition 13. The d- dimensional complementary Grassmann channel Q d is the quantum chan- 
nel defined by the isometry Vg d as G d (ipA' ) = Tr^ o Vg d o ipA' ■ The action of the channel on an 
input audit is given by 

Q\ : VM' h- if = cos 2 ^ 1 ) r tan 2 ( d - fc ) r f d ~ W o x f , (56) 

fe=i ^ ' ' 

where x^P has been introduced in Eq. \21\) . 

For the purpose of proving degradability, W is a harmless unitary matrix independent on k 
so we may just ignore it. But Ork—i)-*(d—k) i s n °t a completely positive map so degradability 
remains to be proven. 

Theorem 14. All Grassmann channels Q d from Eq. are degradable for r £ [0, 7r/4]. 
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r 



FIG. 1: Quantum capacity for Grassmann channels for d — 2,5,10,50 and 100. In the left plot the 
capacity is calculated using the base two logarithm and in the right plot using the base d logarithm. 
The curves follow the order in the legend. 



Remark. Note that the degradability of a qubit erasure channel on the whole interval is recovered 
for d = 2. This corresponds tope [0, 1/2] in Eq. (|23b|) as is valid for a 'standard' qubit erasure 
channel Eq. flU} QJ]. 

Proof. We will prove the theorem by a direct construction of the degrading map. Rewriting 
Eqs. (HJ and ((SSJ) we get 



{d) ~ ©rf = ^PkQdAM (57) 



lA 

k=l k=l 
(1 



T^^V^©?^, (58) 

k=l 

where p k = cos 2 ^ 1 ) r tan 2 ^" 1 ) r(*z\) and Pk = cos 2 ^ 1 ) rt&n 2( - d -V r(j*lj). Recall that for 

the purpose of the degrading map construction we may work with 7^ instead of J^p- We 
assume the existence of the following degrading map 

M : 7 i d) <rf + ± q kX ^ = ±p kX ^ = (59) 

k=1 k=l 

In order for M to be a completely positive map we have to show that < qk and J2k=i 9k = 1 
for re [0,7r/4]. For all d we get from Eq. (1591 the following set of d linear equations 



qipi = Pi 
qiP2 + <?2 = P2 

qiPd + qd=Pd- 



(60) 



The set is easily solvable. The first equation gives us q\ = tan 2 ( d ^ r which we plug into the 
remaining equations. We get 



qk = 



(1 - l) r ( tan2( ^ fc) r ~ tan2(d+fc " 2) r ) (61) 
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for k = 2...d. Since for r G [0, 7r/4] the tangent function is monotonously increasing and 
< tanr < 1 holds as well we may conclude that for k = 1 . . . d all qk are positive. Finally, 
summing the left and right side of the equation set and using Pk = X)fe=i Pk = 1 we find 

that Y?k=i Ik = 1- ■ 

We might proceed to the calculation of the quantum capacity of the Grassmann channels. 
Due to the degradability the quantum capacity formula Eq. @ reduces to the optimized co- 
herent information Eq. (j4]). Furthermore, according to Theorem [7] the Grassmann channels are 
covariant. Therefore, the supremum in Eq. Q is achieved for a maximally mixed input qudit 
and the quantum capacity formula reads 

Q(Q d ) = I cos 2 ^- 1 ) rJ2k( d \ log Jfe( tan 2 ( d - fc > r - tan 2 ^ 1 ) r) . (62) 
k=l ^ ' 

The plots for various d can be found in Fig. Q] where we have plotted the capacities using both 
the base two (on the left) and base d logarithm (on the right). The reason for the presence of 
the base d logarithm is to compare the capacity in a more fair way. 



V. CLASSICAL CAPACITY 



We first present a simple generalization of the characterization theorem derived in [17| (The- 
orem 1). 

Theorem 15 ([17}). The quantum channel T : B(C d ) -> B(C d ') is of the form 

rpi \ Id' - mM(p) 
d' — m 

where M is a positive, linear and trace-preserving map such that there exists a state po where 
mM(po) is a projection of rank m if and only if the a-Renyi minimal output entropy H^ Lin {T) — 
mm p H a (T(p)) is a-independent. 

Remark. The (almost trivial) generalization lies in setting d' ^= d. 

We now show that the blocks from which all the Grassmann channels are composed fulfill 
the required conditions. 

Lemma 16. Every block of the qudit Grassmann channel Qd '■ B(C d ) — > B(C 2 _1 ) has the 
following form: 

I d , - mM d , k (p) 

SdMP) = j, (63) 

d — m 

where d! — (^) is the dimension of the output space. Moreover, for all input pure states po, 
mMd^kipo) is a projection of rank m 

Proof. Consider the fc-th block of the qudit Grassmann channel with input of the form |1)(1|. 
Then the output state will have the following form: 

=7^17 ^ |1 5 2 . . . Sd ><l S2 . . . Sd U (64) 
\k-l) Si 
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where we are using the convention J"\_ 2 s; = k — 1 and normalized the block. Now suppose we 
let m = ( T 1 ) and we know the dimension of the output space for the fc-th block is equal to 
(f) , then d'-m= © - (V) = (£1) • Therefore, 

- (t? - m)g d>fe (|lXl|) = la - I 1 s 2 • • • s <i)(l s 2 . . . s d | (65) 

Si 

= ^|0t 2 ...i d )<0t 2 ...t d |. (66) 

t< 

with the convention now being Yl^—o U — k. The matrix in Eq. (|66l) has rank m = (^fc 1 )- 
Thus with an input of the form po = we have found a matrix Md,k(po) that satisfies the 

condition in Eq. 

Now, in order to generalize this result to arbitrary pure inputs, one can use the SU (d) covari- 
ance of the Grassmann channel presented in Corollary [8] Assume Rd,r to be an r-dimensional 
unitary representation of SU(d). Then due to covariance the following holds: 



Gd,k(P) = Gd,k{Rd,d ° Pa) = Rd,d' ° Gd,k(po)- (67) 
Thus the following chain of equalities hold: 

Rd.d' ° {™>Md,k(po)) = Rd.d' o ((Id' - (d' - m)G d ,k(po)) (68a) 

= I d ' - (d' - m)Rd,d> o Gd,k(Pa) ( 68b ) 

= I d ' - (d' - m)g d M{Rd,d Po) (68c) 

= mM d ,k(Rd,d Po)- (68d) 



The last equality follows from defining mMd,h(Rd,d ° Po) i n this way, since it has rank m since 




FIG. 2: Classical capacity for Grassmann channels for d = 2, 5, 10, 50 and 100. The vertical dashed 
lines at r = 7r/4 label the point where the quantum capacity in Fig. [1] equals zero. In the left plot the 
capacity is calculated using the base two logarithm and in the right plot using the base d logarithm. 
The curves follow the order in the legend. 



it is just equal to a rotated state of rank m. Thus every block, after normalization, has the 
form of Eq. ([531). ■ 



Corollary 17. Grassmann channels have a direct sum form Eq. (|2"Tj) . Hence from Lemma 3 
in [13j it follows that all Grassmann channels have the Holcvo capacity additive. From Eqs. (fTJ) 
and ^ we get 



C{Qd) = sup 



H(A)g d{axAl) -Y,P(x)H(A) gc 



(69) 
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Consequently, exploiting the result from [lj, |3l| for covariant channels, the classical capacity 
is given by 



<d 



C(G d ) = H({p k }) + ]T Pk log - H (A)g d{po 



k=l 

d 



3) 



= log d - cos 2 ^ 1 ) r tan 2 ^" 1 ) r (* _ ^ log k. (70) 



fc=i 



VI. PHYSICAL IMPLICATIONS 

Bosons, fermions and non-inertial observers 

It is often claimed that there is a fundamental difference between the entanglement behavior 
of maximally entangled states built upon fermions and bosons in a relativistic setting [38]. 
Namely, provided that a uniformly accelerated observer has one half of an initially maximally 
entangled state it is found that in the infinite acceleration limit the bosonic entanglement 
disappears while the fermionic entanglement partially persists. Based on this naive approach 
it is concluded that in the infinite limit entangled fermionic states might be useful for various 
quantum-informational protocols where a certain amount of shared entanglement is usually 
needed. 

However, our capacity results suggest something different at least for the purpose of quan- 
tum communication between an inertial and noninertial observer. For the purpose of sending 
quantum messages the ultimate measure of a channel's capability to transmit the information 
is its quantum capacity which has nothing to do with any particular entanglement measure. 
We found that there is no qualitative difference between the behavior of the quantum capa city 
for the qudit Grassmann channel (fermions) and its bosonic equivalent introduced in Ref. [22| 
where we study the bosonic version of transformation Eq. This led to the definition of the 
qudit Unruh channels as the bosonic counterpart of the qudit Grassmann channel. One of the 
resolved problems is the quantum capacity of the qudit Unruh channels which we can compare 
with the Grassmann channels. As a result we find that quantum capacities for both channels 
converge to zero as we approach the infinite acceleration limit. So from the viewpoint of quan- 
tum Shannon theory entangled resources based on bosons or fermions are equally useful. The 
only difference is in the capacity value for a finite acceleration. To fairly compare the capacities 
we rewrite Eq. (|62p as a function of w = tan 2 r. Using cos 2 r = 1/(1 + tan 2 r) we get 

■y / y X ^~ 1 ^ 



QiGd) = d\T+^ ' ^ fc ( «■ 1 log ^ w w > (71a) 



k=l ^ ' 

d^lMVNkr <71b) 



where the details of derivation leading to the second line can be found in Appendix [Bj This 
expression can be directly compared to the one we get for the qudit Unruh channel lAd [22| 

Q(Ud) = -(l- z) d+1 Efc( d+ ^ _1 )log ^pV" 1 . (72) 

The Unruh channel outperforms the Grassmann channel for low-dimensional inputs as we can 
see in Fig. [3l In the infinite acceleration limit z — w — > 1 the ratio of the quantum capacities 
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FIG. 3: Quantum capacity for the Grassmann channel (on the left) and the Unruh channel (on the right) 
taken from j22j] . The parameters w = tan 2 r and z = tanh 2 r are chosen to compare the capacities for 
the same proper acceleration. The capacity of the Unruh channel is in general higher but importantly 
for the infinite acceleration limit (w = z = 1) both capacities are equal zero. The curves follow the 
order in the legend. 



remains finite for all d 

Q(G d ) d\nd 1 L vV,, r^Jd-A, d-k 



- = ss ggg = &r x>- 1 - 2k) { i ) log (73) 



The details of derivation are presented in Appendix [B] We can therefore conclude that the 
capacities asymptotically behave in the same way. We have not found an analytical form for 
the sum in r^ but it is clear that r 2 = In 2 < < r^ +1 for all d. Moreover, numerical simulations 
suggest that r oa = l. 

As a closing comment note that due to the local similarities between the Rindler and 
Schwarzschild spacetime (explicitly spelled out, for example, in [39j ) the similar conclusion 
regarding the quantum capacities also holds in the black hole scenario. 



Grassmann channels and transpose-depolarizing channels 

Let us define the family of qudit transpose-depolarizing channels fl6l ] . 
Definition 18. Let T : B(C d ) — > B(C d ) be a map defined as 

T(a)=ta+(l-t)t/d (74) 

acting on normalized density matrices a. The bar denotes complex conjugation in a given basis 
and the map is a quantum channel in the following interval oft: 

1 <t<^~- (75) 



d-1 ~ ~ d+ 1 

Remark. The transpose-depolarizing channel for t = —l/(d — 1) is known as the d-dimensional 
Werner-Holevo channel [4l| . 

Lemma 19. The complementary channels of the Qd,2 channels presented in Definition are 
the d-dimensional Werner-Holevo channels. 
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Proof. We take the corresponding part of the isometry output Eq. (I19cp for k = 2 and rewrite 
it as 

d d-l 

= c 08 ''" 1 r tanr ^ ^(±)^ |0 . . . l s . . . U . . . 0) A |0 . . . lj . . . 0) c 

i=i 

(76) 

d d-l v ' 



%T3 



Note that = so we can further write (leaving out the irrelevant trigonometric 

functions) 



d d-l 



*ac = E E I m) A ( - A |i> +ft 10 ) c - ( 77 ) 

j=l i>j 

Considering an input pure state \ip} — Yli=i Pi 10-4' i some try of Qd.i reads 

d d-l 

Vs*,2 = E E ( " 5 « I liX^loA' + I liX'IcA') • (78) 



Tracing over A followed by normalizing by Jl/ = y/TJJd — 1) and changing the overall 
sign leads to ( 2 ) Kraus operators of the form 



' K n = \i^j(\m-\m), (79) 

where 1 < i < j < d. These are well known as the Kraus operators for the c?-dimensional 
Werner-Holevo channels [4l|- M 

This result brings us two interesting points. As mentioned earlier, in (To| we study the 
capacity region of the bosonic version of the transformation from Eq. @ known as the qudit 
Unruh channel [22j. One of the results is the characterization of the complementary channels 
of the Unruh channel. Their structure is also block-diagonal as in the present case and the 
first nontrivial complementary block for each d is the transpose-depolarizing channel Eq. (1741) 
for t — 1). This is a peculiar observation and raises a number of questions. First of 

all, why for fermions we get the transpose-depolarizing channel from one end of the allowed 
interval Eq. (fT5j) and for bosons from the other end? Also, does some sort of intermediate 
statistics interpolating between bosons and fermions correspond to the whole interval? One of 
the obvious possibilities are anyons whose appearance is not limited just to the two-dimensional 
world H3- 

The identification of the <i-dimensional Werner-Holevo channel also implies that the Grass- 
mann channels do not belong to the class of Hadamard channels [TJJ, [ll[ - the channels whose 
complementary channel is entanglement-breaking. The reason is that the Werner-Holevo chan- 
nels are known not to be entanglement-breaking. The transpose-depolarizing channels are 
entanglement-breaking for — l/(d 2 — 1) < t [l6j]. Henceforth, at least one of the complementary 
blocks of all Grassmann channels has negative partial transpose. 
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VII. CONCLUSIONS 

We have introduced a new class of quantum channels to the group with computable classical 
and quantum capacities, the Grassmann channels. Such channels are rare in quantum Shannon 
theory since the calculation of the classical and quantum capacities requires an optimization over 
an infinite number of channel uses. The Grassmann channels' isometric extension is physically 
well motivated and stems from the Bogoliubov transformation which preserve the canonical 
anticommutation relations corresponding to the Fermi-Dirac statistics. 

In order to determine the quantum capacity of the Grassmann channels, we have shown that 
these channels are degradable and have explicitly calculated the degrading map. Combining 
this with the result that the Grassmann channels are covariant we were able to give a closed 
form for the quantum capacity. A different technique was used in order to calculate the classical 
capacity of the set of Grassmann channels. We showed that each block, in the block diagonal 
form of the matrix, had a particular form [TtJ which enabled the calculation of their minimum 
output entropy. Exploiting this result allowed us to calculate the classical capacity of the 
Grassmann channels. 

To appreciate the capacity results from the physical point of view we compare the Grassmann 
channels with the Unruh channels studied elsewhere [ill, H3 • They share a close analog to the 
Grassmann channels in the sense that their isometric extension appears formally identical. The 
difference is that for the Unruh channels the isometry is built upon the operators obeying the 
canonical commutation relations (relevant to the Bose-Einstein statistics) inducing a completely 
different class of channels. In particular, the Unruh channels belong to the set of Hadamard 
channels - the channels whose complementary channels are entanglement breaking. The set of 
Grassmann channels do not fall in this class. This is the first example to our knowledge of a 
set of channels that do not belong to the Hadamard class that have a computable and at the 
same time nonzero classical and quantum capacity. 

The main physical consequences also come from the comparison between the fermionic and 
bosonic case. One of the discussed physical motivations for investigating this type of channel is 
that both the fermionic and bosonic Bogoliubov transformation occurs in the study of particle 
production in uniformly accelerated frames. Fermionic entanglement between two parties who 
originally share a maximally entangled state exists to an extent even in the limit of infinite 
acceleration of one of the participants. Contrary to the fermionic case, bosonic entanglement 
vanishes in the infinite acceleration limit. Based on this observation it is believed that there is 
a difference between these types of resources. The result of our work suggests that at least for 
quantum communication purposes there is no difference whatsoever. The Unruh channel does 
demonstrate a greater quantum capacity than the Grassmann channel when the acceleration 
parameter is small, however in the infinite acceleration limit, both quantum capacities tend to 
zero. 

There exists another connection between the Unruh and Grassmann channels. Both channels 
are direct sums of other quantum channels and so are their complementary channels. The first 
non-trivial block of the complementary channel for a given dimension d belongs to the family of 
qudit transpose-depolarizing channels which is a single-parameter family of quantum channels. 
Interestingly, this holds both for the Grassmann channels and Unruh channels. The difference is 
that the complementary block of the Unruh channel corresponds to the transpose-depolarizing 
channel with the parameter at the upper limit of the allowed parameter interval, while the 
first nontrivial complementary block of the Grassmann channel corresponds to the lower limit 
of the allowed interval (the Werner-Holevo channel). This immediately raises the question: 
Why does the fermionic case occupy one end of the interval while the bosonic case occupy the 
other? Perhaps and even more interestingly, could there be intermediate statistics model (like 
anyons, for example) that would explain intermediate values of the interval? Another direction 
in which future research could be done is to consider the case of coupled fermions with a 
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more active role of the spin variable and one can ask how the additional spin variable alters 
the Grassmann channels and whether the classical and quantum capacity is still calculable. 
Finally, the implications of the results here obtained to quantum information protocols inspired 
by Cooper pairs in solid state or atomic physics scenarios deserve an independent detailed study. 
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Appendix A: Geometric picture of the sl(d, C) Lie algebra representations 

In the s\(d, C) case there may be several mutually commuting operators. A maximal linearly- 
independent commuting set of operators of a (semi-simple) Lie algebra is called a Cartan sub- 
algebra. Once a Cartan subalgebra has been chosen we can use the common eigenvectors to 
label the basis vectors of an irreducible representation. 

Definition 20. An r-tuple a — (a%, . . . ,a r ) of complex numbers is called a root if: (i) not all 

the cti are zero, (ii) there is an element E a ofs\(d,C) such that 

[H i ,E a ]=a i E a . (Al) 

The set {Hi, E a } is called the Cartan- Weyl basis [4(J. The Cartan subalgebra of sl(d, C) has 
dimension r = d — 1; we say that the rank of the Lie algebra is r. We write H = (Hi, . . . , H r ) 
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for the Cartan subalgebra generated by the elements {Hi, . . . ,H r } of the Lie algebra; these 
elements are assumed to be independent. Among all the roots there is a class of special roots 
called simple roots. 

Definition 21. A root is called a simple root if it cannot be written as a linear combination of 
other positive roots. 

Definition 22. If p (where p : s\(d, C) — > GL(V) for some V) is a representation of s\(d, C) 
then a r -tuple p — (pi, . . . , p r ) of complex numbers is a weight for p if there is a nonzero vector 
ip G V such that ?p is an eigenvector of each Hi with eigenvalue pi . 

If p is a weight and ip a weight vector for p and a is a root then 

p{H t )p{E)iP = (p t + a t )p{E)i>. (A2) 

In short, E changes all the eigenvalues of the Cartan operators and it creates a new weight 
vector (or kills the weight vector). The root is a vector in the weight space that points in 
the direction in which the weights are changing. The E operators are called shift or raising 
and lowering operators. Roughly speaking, the positive roots correspond to raising operators 
while the negative roots to lowering operators. We classify the irreducible representations by 
the highest possible value of the weight. For general semi-simple Lie algebras we do exactly the 
same thing once we have a suitable order on the weights in order to define the right notion of 
highest weight. 

A special example of the Cartan- Weyl basis is the Chevalley-Serre basis [40]. Two aspects 
make this basis special, (i) The step operators are associated to simple roots and (ii) the 
normalization is chosen such that the roots are integers. Unless explicitly stated we work in 
this basis due to its accessible geometric interpretation. 

To give the operators a geometric interpretation we will work with a specific matrix repre- 
sentation of the sl(d, C) algebra. Let us define Eij, where 1 < i ^ j < d, as the matrix having 
one where the i-th row and the j-th column intersect and the rest are zeros. Furthermore we 
define a diagonal matrix Hi in which the i-th diagonal entry is 1, the (i + l)-th diagonal entry 
is —1 and the rest are zeros. If we assume j = i + 1 the following set {Hi, Eij, EjA forms the 
Chevalley-Serre basis for sl(d, C). More explicitly, for d — 2 we get 

Hi = (I (A3a) 

E12 = (° Q (A3b) 

E\ 2 = (J o) ■ ( A3c ) 

Therefore in this basis we have 



[H i ,E ij ] = 2E ij (A4a) 
[H i ,Ej j ] = -2E} j . (A4b) 

The simple roots are elements of a vector space dual to the one spanned by elements of the 
Cartan subalgebra. 

This structure opens the door to an insightful geometric picture in terms of the so-called root 
space diagram. The dual space will be called the space of roots. For the s\(d, C) Lie algebra the 
space is (d — l)-dimensional and the simple root vectors defined as an r-tuple of simple roots 
form a non-orthogonal basis. Let us call the basis spanning the space of roots the root basis. 
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n'"=(o,-D ^'=(0,-1,-1) 



a) b) 

FIG. 4: (a) We illustrate the spaces of roots for sl(2, C) and si (3, C). The basis vectors are simple roots 
and are indicated by a?;. For the weights (vertices) we explicitly write down their coordinates in this 
basis. We can read off the generators of the Cartan subalgebra from the weights, (b) It is sometimes 
helpful to introduce an overcomplete basis. The position of 03 reflects the relation H3 = Hi + H2 for 
Hz in Eq. (|A6[) and it is not a simple root. 



We are aware of the overuse of the expressions root and root vectors. The terminology is not 
stabilized and differ in various textbooks. Also, the root space diagrams are sometimes called 
weight diagrams. It can be shown that each consecutive simple root vectors subtend the angle 
27r/3. We rewrite Eq. (|A"Tj) as 

[H ii E ii ] = {j l ^-^)E ij , (A5) 

(i) 

where fj,) are called fundamental weights. The word fundamental reflects the fact that we are 
dealing only with simple roots. We explicitly rewrite Eq. (| A2|) as the spectral decomposition 



where j labels the j-th component of a d-component spinor | ip} . The important role played by 
the (fundamental) weights is that they are coordinates of the eigenvectors in the space of roots. 
Fig. [4^, illustrates the situation for d = 2 and d = 3. Hence for each fundamental representation 
of s\(d, C) there is d points each representing an eigenvector i[>j- 

What about the role played by the shift operators? They have a precise geometric interpre- 
tation as well. All points of the sl(d, C) fundamental representations are interconnected. The 
operator responsible for a transition from site to \4>j) is the operator E^ or Ejj for the 



opposite direction. It also follows from Eqs. (|A4[) that each segment connecting two neighboring 
spinors has length two. 

The fundamental representations of sl(d,C) algebra contain r — d — 1 independent sl(2, C) 
subalgebras each satisfying Eqs. (|A3|) . However, since the root space diagram is a complete 
graph there are in total ( 2 ) linearly dependent si (2, C) subalgebras corresponding to the num- 
ber of edges. The diagonal generators of the 'additional' sl(2,C) subalgebras are constructed 
similarly to the H^s above. The only difference is that 1 and —1 on the diagonal are separated 
by one or more zeros. As an example id = 3), the remaining diagonal generator is 



(A6) 
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Therefore, the rank of the weight vectors equals three and it correspond to introducing an 
overcomplete root basis, see Fig. Hp. Note that H% does not correspond to a simple root. 
Indeed, the axis a.3 in Fig. 0b can be obtained by a linear combination of a\ and ai which are 
both positive root vectors. 

Completely antisymmetric representations of the sl(d, C) algebra 

For the purpose of this article we are interested in particular higher-dimensional representa- 
tions of si (<i, C) - the completely antisymmetric representations. The lowest-dimensional anti- 
symmetric representation is formed as the dual of the fundamental representation 

= £ A~.i^. 2 ... fe) (A7) 

where e ,1 '" 3d is a completely antisymmetric tensor and 1 < < d. The dual 

representation (of the fundamental representation only!) coincides with its complex conju- 
gate representation. We get the dual representation from the fundamental representation by 
G s \(d,c) ^ — G s \(d.c) where G s i(d,c) are all generators of sl(d,C). The eigenvalues remain the 
same and the weight vectors just change the sign. Hence, the root space diagrams are the same 
and they are just, vaguely speaking, pointing in the opposite direction. The interpretation of 
the edges and points follows the fundamental case. The difference lies in the fact that the 
antisymmetric representations of sl(d, C) are carried by c£-component antisymmetrized spinors. 

Higher-dimensional completely antisymmetric representations of the sl(d, C) algebra are 
formed in an intuitive way. First of all, the dimension of the spaces of roots remains the 
same as well as the number of algebra generators. The generators clearly satisfy the same com- 
mutation relations since it is just a different representation of the same sl(cZ, C) algebra. The 
root diagram 'is grown' in the direction of roots but this process cannot go on forever. The 
spinors carrying the higher-dimensional representation are completely antisymmetrized and so 
there are only (^) states in the fc-th completely antisymmetric representation of s\(d, C) where 
1 < k < d (formally, we should also include the trivial representation, k = 0). Fig. [5] illustrates 




a) b) 

FIG. 5: (a) The left plot is the dual representation to the fundamental representation of si (3, C) 
where ip l are completely antisymmetric spinors from Eq. (|A7|) . The middle plot is the second- lowest 
(k = 2) completely antisymmetric representation of si (3, C). Due to the antisymmetrization procedure 
(indicated by the square brackets) the vertices do not correspond to any state, (b) The second lowest 
(k = 2) completely antisymmetric representation of si (4, C). The antisymmetrization procedure leaves 
us with a six-dimensional space (the inner octahedron) with the spanning basis indicated. The spinors 
are written in the Fock basis following Definition [1] The sign attached to an edge indicates the sign 
induced by the (fermionic) sl(2,C) shift operator when applied on spinors the edge connects. 



the growth of these representations for d = 3 and d = 4. The connection to the fermionic states 
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brought in Definition [I] is straightforward: The state \F} = \n\ . . .rid) for which X^=i n>i = k 
holds is an antisymmetric spinor carrying the A:-th completely antisymmetric representation of 
sl(d,C). 

We find the matrix form of the generators of all higher-dimensional completely antisym- 
metric s\(d, C) algebra representations. Because of the complete antisymmetry condition there 
are only segments connecting two neighboring points. Therefore only the lowest-dimensional 
representation of the sl(2,C) subalgebra appears in the construction. However, because for a 
given edge there might be more segments parallel to it, the sl(d, C) subalgebra generators are 
formed by a direct sum of the sl(2, C) subalgebras. The reason for a direct sum is that they, by 
construction, act on mutually orthogonal subspaces. Finally, for every d there is only r = d — 1 
linearly independent directions (or, said otherwise, only r — d — 1 independent sets of parallel 
lines) so the number of generators equals the number of independent sl(2, C) subalgebras and 
they manifestly satisfy the commutation relations for sl(d, C). Note that it does NOT mean 
that the completely antisymmetric representations of sl(rf, C) are direct sum representations. 
As a matter of fact, they are irreducible. The direct sum subalgebras we have created do not 
themselves span mutually orthogonal subspaces. For illustration see Fig. [5}) where any pair of 
parallel segments 'share' spinors with some other pair of parallel segments. 

The only ambiguity lies in the sign of the shift operators. We can see from Eq. (|A3[) that 
switching their sign does not spoil the commutation relations. Let us illustrate the ambi- 
guity on an example from Fig. [SJj. The transition from state 1 001 1) to |010l) is provided 
by a\a,3 |001l) = |010l) meanwhile to get from |001l) to |0110) we acquire a minus sign 
a\ci4 |001l) = — |0110). Not accidentally, the operators ajcij,aj<ij are the fermionic representa- 
tion of the shift operators of the si (2, C) algebra po| . They play an important role in the proof 
of Theorem [7] where the fermionic representation is properly introduced. 

So far we talked about a direct sum of several si (2, C) subalgebras but for the sake of proof 
of Theorem [7] we need to specify how many summands there actually is. This transforms 
into the question how many different sl(2,C) subalgebras corresponding to a chosen direction 
exist. Let the direction be chosen by the step operators {a\ai,a\a,2\. Then for a given d and 
k there is (fl?) parallel segments in this particular direction for the fc-th lowest completely 
antisymmetric representation of sl(d, C). We get this number by a combinatorial argument: We 
have a spinor with d positions where the first two slots are occupied. For the lowest-dimensional 
representation (A; = 1) the rest is occupied by zeros and therefore for the k-th lowest dimensional 
representation there is (k — 1) ones to occupy the remaining (d — 2) free spaces. By a simple 
permutation argument we can see that this holds for any of the (d — 1) linearly independent 
directions. For the example in Fig. [5}d we indeed get two parallel segments for each of the three 
independent directions. 



31 



Appendix B: The capacity ratio derivation 

On can show that Eq. (|71al) simplifies in the following way: 



QM = ^TT^^" ~ w " 1)fc U) logfc (B1) 

k— 1 x / 



i^p+r (E ( d - (/_ z ) lo § ( rf - - E © l0 § fc ) ( B2 ) 
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Now consider the Unruh channel quantum capacity from Eq. (|72p . 



^, 7( , (l-z) d+1 ^ , , fd+k- 1\ , d + fc-1 , . 
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Now as z — > 1 the dominant terms in this summation are the terms such that k 3> 1. Thus in 
the limit z->lwe can approximate the log using the following: 

, d + k - 1 , / d - 1 \ 1 d - 1 , . 



Let us recall that log is the logarithm base d and In denotes the natural logarithm. Eq. (|B6I) 
becomes: 

d lnd z(l-z) d 1 ' 

am a z 

Finally, for w = z — > 1 the ratio reads 

QiGd) ,. dlnd z ^ i. j_ 2 _ 2 t. (d — 1\ , d— k 

Urn ~!, i = lim 7 tt-tt ; r^r > z k (l + z-\ z d 2 lk )[ log 
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The following lemma makes sure that the approximation in Eq. (|B7[) is valid. 
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Lemma 23. The difference \Q(Ud) — Q'(Ud)\ approaches zero at a rate O ((1 — z) 2 ) for z — > 1. 
Proof. We can split Eq. (|B6[) to approximate the logarithm, 
rd-i 
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For k > d we can Taylor expand the logarithm, 
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Since this is a converging alternating series, we can bound the following quantity for k > d, 
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Consider the following approximation of the quantum capacity (rewritten Eq. (|B8p ) 
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